Field-induced phase transition in the periodic Anderson model 
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Abstract 

We investigate the effect of magnetic fields on a Kondo insulator by using the periodic Anderson model. The 
analysis by dynamical mean field theory combined with quantum Monte Carlo simulations reveals that the magnetic 
field drives the Kondo insulator to a transverse antiferromagnetic insulator at low temperatures. We calculate the 
staggered spin susceptibility and find its divergence signaling the antiferromagnetic instability. Further investigation 
of the spin correlation functions and the magnetization process clarifies how the magnetic field suppresses the 
Kondo singlet formation and induces the transverse antiferromagnetic ordering. 
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Heavy-fermion systems with various ground states 
have attracted continued interest. One of the inter- 
esting examples is the Kondo insulator (KI), where 
a small charge gap renormalized by the Coulomb in- 
teraction appears at low temperatures [1]. ft is known 
that this insulating phase gets unstable upon introduc- 
ing the magnetic field, the pressure, etc. Experiments 
on some KI in high magnetic fields indicate closure of 
the Kondo- insulating gap [2], exemplifying a transition 
from the KI to a correlated metal[3,4]. If the KI is in 
the proximity of magnetic instability, the local singlet 
formation gets weak, and an applied field may possi- 
bly trigger a phase transition to the antiferromagnetic 
(AF) ordered state before it becomes a metal. 

The periodic Anderson model (PAM) at half filling 
may be a simplified model to describe the Kondo in- 
sulating phase and the AF phase. The magnetic in- 
stability of the PAM has been investigated by a va- 
riety of methods[5,6,7,8]. Recently, the magnetic-field 
effects on the two-dimensional Kondo lattice model 
have been studied on the basis of the mean field theory 
and quantum Monte Carlo (QMC) simulations. It has 
been found that the magnetic field induces a second- 
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order phase transition from the paramagnetic to the 
AF ground state[9,10]. 

In this paper, we investigate field-induced AF phase 
transitions in the PAM by using dynamical mean field 
theory (DMFT) [II] combined with QMC simulations 
[12] at finite temperatures. The model we study here 
is the PAM with the Zeeman splitting, 
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where all the symbols have their usual meaning. We 
consider the hyper-cubic lattice with a bipartite prop- 
erty in infinite dimensions. The bare density of states 
for conduction (c) electrons is Gaussian with the width 



f:po(e)=expt-(e/f) 4 ]A 



We choose t* — 1 as 



the energy unit. The magnetic field B applied along 
the z direction is coupled to both of c and / electrons. 
We set gfiB = 1 for both c and / electrons. 

DMFT is a powerful framework to study strongly 
correlated electron systems. In DMFT, the original lat- 
tice model is mapped onto an effective impurity model 
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Fig. 1. The inverse of the staggered spin susceptibili- 
ties. The dashed lines are fitting functions with the form 
XxxCQr 1 tx B — B c - 
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with the self-consistently determined medium, which is 
solved by means of QMC simulations in our approach. 
We use the typical parameters U = 2.0 and V = 0.6 
in the following calculation. For these parameters, the 
ground state at zero field is a paramagnetic KI, as 
shown by Jarrell et al. [7]. 

We first show that the magnetic field induces a KI- 
AF phase transition. We calculate the staggered spin 
susceptibility Xxx(<l = Q), with Q = [tv,tv,. . .], where 
the suffix x denotes the direction perpendicular to the 
field. In Fig. 1, we plot the inverse of the total suscep- 
tibility Xxx (Q) an d the susceptibility for / electrons 
Xxx(Q) as a function of the magnetic field for differ- 
ent temperatures. At temperatures T = 1/24 and T = 
1/28, the staggered susceptibilities take a maximum at 
the magnetic field B ~ 0.28. At the lowest tempera- 
ture T = 1/32 the susceptibilities diverge at the field 
0.25 < B < 0.34. The magnetic field triggers a phase 
transition from the paramagnetic KI to AF ordered 
state. We determine two critical values of the field by 
extrapolating the inverse susceptibility to zero with the 
form Xxx(Q)' 1 oc B — B c . 

We further calculate the magnetization process and 
the field dependence of the spin correlation functions 
at T — 1/30, which is slightly higher than AF tran- 
sition temperature. The magnetization M, defined as 
M — ([rj| +nt] — [n\ +n{]), is shown in the top panel of 
Fig. 2. We can see the crossover behavior from the KI 
with the spin gap to the paramagnetic metal around 
B ~ 0.28. Around this field, we also see the change in 
the character of the spin correlation functions. In the 
middle of Fig. 2, we show the variance of the /-moment 
in the z direction ([: S% :] 2 , where :: denotes the nor- 
mal order, and in the x direction ([Si] 2 }. When the 
magnetic field larger than B ~ 0.28 is applied, spin 
fluctuations in the z direction are suppressed and the 
moment in the x direction get large, since the Kondo 
singlet formation is suppressed by a magnetic field. 
The consistent behavior is seen in the spin correlation 
functions between the c and / electrons, (: S° :: S[ :) 



Fig. 2. The magnetization process (top), the field dependence 
of the variance of the spin moment (middle) and the spin cor- 
relation functions between the c and the / electrons (bottom). 

and (SxSl), plotted in the bottom of Fig. 2. AF cor- 
relations between the c and / electrons are suppressed 
by a magnetic field. We can see that the correlations 
in the z direction are more suppressed than those in 
the x direction. The above behavior explains why the 
transverse rather than longitudinal antiferromagnetic 
ordering is favored in finite magnetic fields. 
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